Renormalization-group exponents for superconducting phases in two-leg ladders 
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In previous studies, we proposed a scaling ansatz for electron-electron interactions under renor- 
malization group transformation. With the inclusion of phonon-mediated interactions, we show that 
the scaling ansatz, characterized by the divergent logarithmic length Id and a set of renormalization- 
group exponents, also works rather well. The superconducting phases in a doped two-leg ladder are 
studied and classified by these renormalization-group exponents as demonstration. Finally, non- 
trivial constraints among the exponents are derived and explained. 
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Even though superconductivity has been observed [T] 
for one hundred years, the mechanism turning electrons 
into pairs remains an open question. In conventional su- 
perconductors, phonons mediate effective attractions [5] 
between electron in the low-energy limit and lead to 
Cooper pair formation near the Fermi surface. On the 
other hand, for unconventional superconductors^ like 
cuprates, phonon-mediated interactions seem to play a 
secondary role while the electronic correlations are be- 
lieved to reigns. It is speculated that the spin fluctua- 
tions is the key to the pairing mechanism (and perhaps 
the pairing symmetry as well). However, it hasn't been 
fully understood how Coulomb repulsion eventually glues 
electrons into pairs. 

The discovery of iron-based superconductors keeps the 
puzzling charm going. Collecting from experimental ob- 
servations, it is believed that electronic correlations in 
these materials are important but much weaker than 
those in cuprates. Though the renormalization-group 
(RG) analysis 4-6 for electron-electron interactions de- 
livers the correct pairing symmetry, the isotope effects 
measured in laboratories show conflicting results. In ad- 
dition, the iron-based superconductors are distinct from 
the cuprates due to the presence of multiple active bands. 

Inspired by the anomalous isotope effects in iron-based 
superconductors, it is helpful to investigate the competi- 
tions between electron-electron[71[S| and electron-phonon 
interactions [HHH] by RG analysis. The major difficulty 
lies in the enormous couplings for all allowed interactions. 
Besides, as the number of couplings grow, reading out the 
desired messages from RG flows can be challenging as 
well. Recently, we found a scaling ansatz [H], character- 
ized by a set of RG exponents, for electron-electron inter- 
actions in many correlated systems. These RG exponents 
build a clear hierarchy of relevant couplings and serve 
as an unambiguous indicator for the ground-state insta- 
bilities. We are curious whether similar scaling ansatz 
exists even when the retarded interactions mediated by 
phonons are included. 

In this Rapid Communications, we show that the scal- 
ing ansatz works in the presence of both electron-electron 
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FIG. 1: (a) Band structure for a two-leg ladder at filling < 
n < 1. The low-energy excitations are described by four pairs 
of chiral fields ipma, tpLia with i = 1, 2 and a =t,4^- (b) Phase 
diagram for a doped two-leg ladder. In the filling regime 
0.5 < n < 0.54, the interband Cooper scattering becomes 
irrelevant and the superconducting phase decomposes into the 
2-channel Luttinger liquid. 

and electron-phonon interactions and the patterns of 
the extracted RG exponents describe the ground states 
and also the quantum phase transitions between them. 
To make the discussions concrete, we study the exem- 
plar two-leg ladder shown in Fig. 1. The low-energy 
excitations are described by four pairs of chiral fields 
^Ricn V'tiaj where R/ L denotes the chirality and i = 1, 2, 
a =t, i- are the band and spin indices. The generic inter- 
actions can be classified into two categories: forward and 
Cooper, described by the following effective action 

121, 

where and Cij denote the forward and Cooper scatter- 
ing between the i—th and j— th bands. To avoid double 
counting, the diagonal parts of the forward scattering are 
set to zero, fa = 0. The superscripts I, s stand for large 
and small momentum transfer between chiral fields of the 
same spin index. 
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The ground state of the doped two-leg ladderjTH [TH 
[22] depends on the fiUing n, i.e. average number of 
electrons per lattice site. For 0.54 < n < 1, electrons 
form pairs and the ground state is an unconventional su- 
perconductor (with power- law correlations). The wave 
functions in different bands reveal opposite signs, resem- 
bling the pairing symmetry in iron-based superconduc- 
tors. That is to say, the pairing wave functions at fc^ = 
(band 2) and ky — tt (band 1) are opposite in signs and 
thus referred as "d-wave" unconventional superconduc- 
tor (d-USC). In the filling interval, 0.5 < n< 0.54, the 
pairing instability disappears and the ground state is a 
2-channel Luttinger liquid. Bellow the filling n < 0.5, 
only one band remains active and the low-energy physics 
is well described by the single-channel Luttinger liquid. 
Despite the simplicity of the phase diagram, it is worth 
mentioning that RG exponents for the relevant couplings 
help to clear up some confusions in the literature. 

Electron-phonon interactions are different in nature 
and retardation effects must be taken care of. To in- 
clude both electron-electron and electron-phonon inter- 
actions, it requires turing the couplings gi into frequency- 
dependent functions gi{uj), making the RG analysis chal- 
lengingly difficult. The old trick in the original Bardeen- 
Cooper-Schrieffer (BCS) theory comes to rescue - we 
assume the phonon-mediated interactions only occur 
within a thin shell, roughly the order of Debye frequency 
ljd, near the Fermi surface. Within this approxima- 
tion, the frequency-dependent couplings take the follow- 
ing form, 

gd^^) ^ 9i + ^{^D - i^)9i, (2) 

where Qi and tji are the (instantaneous) electron-electron 
and the (retarded) phonon-mediated interactions. The 
step-like shape is invariant under RG transformation. Let 
1 and O represent the shapes of constant and step respec- 
tively. It is clear that the algebra of the shapes is closed: 
1 X 1 = 1, 1 X e = e X 1 = e and e X e = e. As a con- 
sequence, the renormalized couplings maintain the same 
step-like shapes and make the approximate frequency de- 
pendence self-consistent. 

Following standard derivations, the RG equations for 
the retarded interactions are 

d,c% = - 4i^4 - 44i^ + 2c^/^ - 442^ + 2i- 5^, 
- - 2(i^)2 - 2(4)2 - 44/^ + 2i- 

-44£f^+24cf^. (3) 

The other set of retarded interactions cf, , does not 
renormalized at all within the approximation. The renor- 
malized couplings gi{l) and gi{l) are solved numerically. 
It was found before that gi{l) is well captured by a scal- 
ing ansatz. It is rather remarkable that the numerical 
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FIG. 2: Phase diagram for a doped two-leg ladder with both 
electron-electron and electron-phonon interactions. The crit- 
ical filling He « 0.57 separates two types of superconduct- 
ing instabilities driven by electron-electron interactions and 
phonon-mediated ones. 

solutions for gi also follow a similar scaling ansatz, 



where Id is the divergent (logarithmic) length, jg. is the 
RG exponent for the retarded coupling gi and Gi is some 
non- universal order-one constant. We shall elaborate the 
importance of the divergent length Id and the RG ex- 
ponents in later paragraphs. It is worth emphasizing 
that all renormalized couplings remain in the perturba- 
tive regime even though the ansatz mysteriously contains 
the divergent length scale Id- 

For numerical analysis, we choose the ratio gi/gi = 
—0.1, where the minus sign indicates the attractive na- 
ture of the phonon-mediated interactions. As shown in 
Fig. 2, an interesting critical filling Uc ~ 0.57 separates 
two types of superconductors driven by different pairing 
mechanisms. In the regime ric < n < 1, the instanta- 
neous couplings gi{l) grow faster than the retarded ones 
gi{l) and the divergent length Id is dictated by the bare 
values of the electron-electron interactions. Compared 
with the phase diagram in Fig. 1, the superconducting 
ground state is qualitatively the same with unconven- 
tional pairing symmetry between different bands. In the 
regime < n < Uc, the RG flows are different. Despite of 
the smaller bare values, the retarded couplings gi{l) reign 
and flow toward strong coupling first. The pairing mecha- 
nism is the well-known BGS theory and the isotope effect 
is expected to be significant. This regime is further di- 
vided into two parts because there is only one active band 
at lower fillings. In the two-band regime, 0.5 < n < ric, 
Van Hove singularity develops in the almost empty band 
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these superconducting phases. In the fiUing range 0.57 < 
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FIG. 3: RG exponents of intraband Cooper scattering in dif- 
ferent superconducting phases (a) for the anti bonding band 
(b) for the bonding band. The pattern of the RG exponents 
provides a simple yet unique indicator to classify different 
ground states. 



and the phonon-mediated Cooper scattering (cn) drives 
the ground state superconducting. Compare with the 
phase diagram in Fig. 1, inclusion of the electron-phonon 
interactions causes the instability in the 2-channel Lut- 
tinger liquid and only one channel (in the majority band) 
survives. Further reduction in electron filling to the 
regime n < 0.5, only the bottom band remains active. 
This is where our common intuitions work - the repulsive 
electron-electron interactions (though with larger bare 
values) renormalize to zero while the attractive interac- 
tions mediated by phonons reign, rendering the ground 
state to the conventional BCS superconductor. 

The critical density Uc not only separates two different 
pairing mechanisms, but also marks different trends of 
the divergent length Id at different fillings. In the regime 
< n < ric, Id remains more or less constant when the 
filling varies. The constancy partially comes from the 
smooth evolution of Fermi velocities in this regime and 
is partially attributed to the steady changes of the inter- 
band Cooper scattering. On the other hand, the diver- 
gent length Id shows sensitive dependence in the regime 
< n < Uc- This is expected because the (mean-field) 
transition temperature is roughly Tc ~ e"''' by standard 
scaling argument. Compared with BCS theory, Id is in- 
verse proportional to the density of states and thus pro- 
portional to the Fermi velocity (in quasi-one dimension) . 
The almost-linear dependence of Id in Fig. 2 just reveals 
the underlying velocity evolution. 

Now we come to the RG exponents (shown in Fig. 3) in 



n < 1, the RG exponents: 7c 
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among all and the pairing mechanism mainly arises from 
the electron-electron interactions. The RG exponent 
is around 0.6 near half filling n = 1 and remains almost 
constant (but with slight decrease) upon filling reduction. 
On the other hand, 7522 starts around 0.6 near n = I but 
increases to almost unity when approaching n « 0.57. 
The enhanced retarded Cooper scattering in the bonding 
band 622 correlates with the larger spin gap from the 
electronic correlations. Although the phonon-mediated 
interactions are also relevant (positive RG exponents), 
the exponents 7£ii,7£22 the retarded interactions are 
smaller than one and only play a secondary role in this 
regime. In short, one can say that the pairing mechanism 
in this regime is mainly attributed to electron-electron 
interactions. However, the observed physical quantities 
are "dressed up" by retarded couplings and may depend 
on the electron-phonon interactions significantly. 

With fining 0.5 < ?i < 0.57, BCS instability occurs in 
the nearly empty band due to the Van Hove singularity. 
Note that this regime is almost twice larger than the 2- 
channel Luttinger hquid phase (0.5 < n < 0.54) when 
the phonon-mediated interactions are ignored. As shown 
in Fig. 3, all RG exponents in this regime are zero ex- 
cept jcii = 1- The pattern of the RG exponents makes 
the analysis rather straightforward. Electrons form pairs 
driven by the effective attraction in the (almost empty) 
antibonding band and other couplings remain small even 
at the cutoff length scale. The pairing mechanism is fur- 
ther supported by the linear dependence of the divergent 
length Id (X {n — ric) seen in Fig. 2. It is computed that 
the density of states in this regime is mainly ascribed to 
the antibonding band and the emergent Van Hove singu- 
larity takes the form l/{n — nc). Following BCS analysis, 
the critical temperature for the superconducting phase is 
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On the other hand, scaling 
argument gives Tc ^ e"'"*. Comparing both approaches 
together, the linear dependence of the divergent length 
Id is explained. 

Entering the one-band regime, < n < 0.5, the pattern 
of the exponents are also quite simple. As shown in Fig. 
3, all RG exponents in this regime are zero except 7522 = 
1. The pairing instability now switches to the bonding 
band and the effective action of the ground state is well 
captured by BCS theory. The overall linear trend of the 
divergent length Id (shown in Fig. 2) follows the inverse 
of the density of the states as explained in the previous 
paragraph. 

The pattern of RG exponents build up the hierarchy of 
relevant couplings and serves as an unambiguous indica- 
tor of the ground state. Furthermore, the pattern change 
is quite sharp when going through a quantum phase tran- 
sition. The usage of these exponents are not limited to 
the simple two-leg system. We have applied the classifi- 
cation scheme to more general correlated systems in two 
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dimensions and the primitive results suggest the RG ex- 
ponents can also be extracted numerically. In fact, the 
proposed scaling ansatz has a solid root from mathemat- 
ical aspect. For typical divergent flows studied here, a 
so-called ■!/'-series[53] is shown to exist with a single di- 
vergent length scale Id- The RG exponent represents the 
leading order contribution of the infinite series expansion. 
Apparently, in most parameter regimes, the leading order 
term seems to grab the essential information in the RG 
flows. Due to the flexibility and generality of the ?/;-series 
expansion, it is not yet clear whether the leading-order 
term can always be made dominant. At presence, the 
dominance of the leading order term is demonstrated in 
numerics without rigorous proof. 

The dominant scaling ansatz also gives rise to non- 
trivial constraints between the RG exponents. For in- 
stance, plugging scaling ansatz for phonon-mediated in- 
traband Cooper scattering cu in Eq. ([s]), it would imply 
either both cu and cu are irrelevant at the same time, 
or the exponents must satisfy the constraint, 75.. -I- 1 = 
max{27g. . , 7c.; +75;;}. That is to say, the RG exponents 
for the intraband Cooper scattering are constrained, 

7£ii=7cH=0, or max{7g..,7c,J = 1. (5) 

These constraints are apparent in the extracted RG ex- 
ponents shown in Fig. 3. Other constraints on the RG 
exponents can be derived in similar fashion and agree 
with our numerical results at hand. Finally, we like to 
conclude that the scaling ansatz characterized by the di- 
vergent length Id and the RG exponents 7i also works for 
phonon-mediated interactions. The classification scheme 
proposed here is simple but powerful without ambiguity. 
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